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Abstract 

It is now commonly accepted that cosmic 7 -ray bursts (GRBs) 
are of cosmological origin. This conclusion is based on the statisti- 
cal analysis of GRBs and the measurements of line redshifts in GRB 
optical afterglows, i.e., in the so-called long GRBs. In this review, 
the models of radiation and models of GRB sources are considered. 
In most of these models, if not in all of them, the isotropic radiation 
cannot provide the energy release necessary for the appearance of a 
cosmological GRB. No correlation is noted between the redshift, the 
GRB-spectrum shape, and the total detected energy. The compari- 
son between data obtained in the Soviet experiment KONUS and the 
American experiment BATSE shows that they substantially differ in 
statistical properties and the detection of hard x-ray lines. The in- 
vestigation of hard gamma (0.110 GeV) afterglows, the measurement 
of prompt optical spectra during the GRB detection, and the further 
investigation of hard x-ray lines is of obvious importance for gaining 
insight into the GRB origin. Observations of two bright optical GRB 
afterglows point to the fact that an initially bright optical flare is di- 
rectly related to the GRB itself, and the subsequent weak and much 
more continuous optical radiation is of a different nature. The results 
of observations of optical GRB afterglows are discussed. They point to 
the fact that the GRBs originate in distant galaxies with a high matter 
density, where intense star formation takes place. The interaction of 
the cosmological GRB radiation with a dense surrounding molecular 
cloud results in the appearance of long-duration (up to 10 years) weak 
optical afterglows associated with the heating and reradiation of gas. 
Results of 2D numerical simulation of the heating and reradiation of 
gas in various variants of the relative disposition of GRB and molecu- 
lar clouds are presented. In conclusion, the possible relation between 
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the so-called short GRBs and recurrent sources of soft 7 rays in our 
Galaxy, the so-called soft gamma repeaters, is discussed. 



1 Introduction 

Nowadays, it is commonly accepted that cosmic 7 -ray bursts (GRBs), the 
discovery of which was reported in 1973 [55], are of cosmological origin. The 
first cosmological model based on explosions in active galaxy nuclei (AGN) 
was proposed in [84]. GRB formation near a collapsing object due to neutrino 
antineutrino annihilation was investigated in [12]. Previously, the model of 
GRB formation during collapses and explosions of supernovae was considered 
in [23], where the following chain of reactions resulting in gamma- photon 
production was considered: 



In brackets, wc give photon energies. The reactions with heavy nuclei were 
also considered in [23]: 



The starquake, the subsequent explosion, and the ejection of matter from 
a nonequilibrium layer in the neutron-star crust discovered in [22] are accom- 
panied by the gamma rays induced by the fission of superheavy nuclei. This 
scenario was proposed in [23] as an alternative model of the galactic GRB 
origin. The total number of various galactic models exceeds one hundred. 
Even now, after the discovery of rcdshifts in optical afterglows, which became 
possible due to the x -ray observations of GRB with the Beppo-SAX satellite 
and their subsequent optical identification, the galactic models are of more 
than simply historical interest because the origin of short GRBs with the 
duration of less than 2 s remains uncertain. In [12], it was obtained that the 
efficiency of transformation of neutrino-flux energy ~ 6 ■ 10^^ erg into the 
energy of the X - and 7 -ray burst amounts to a fraction a ~ 6 ■ 10^ with a 
total GRB energy release Wx,-y ~ 3 • 10^^ erg. The 3D numerical simulation 
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of two colliding neutron stars [86] and a hot torus around a black hole [87] 
showed higher efficiency oi X - and 7 -ray production, reaching 0.5% in the 
first case and 1% in the second case. In part, the distinction with [12] can be 
associated with the geometry more preferential for the neutrino-flux out- flow 
when the annihilation rate increases in comparison with that estimated in 
the spherical geometry [12]. Nevertheless, even in such an optimistic variant, 
the GRB formation with a total energy yield in the X - and 7 -ray regions not 
exceeding 5 • lO^'^ erg is probable. This energy is insufficient for explaining 
the energy resource of many GRBs because only the direct optical GRB 
radiation can attain 10^^ erg, and the isotropic flux in the gamma region 
reaches 2.3 • 10^^ erg for GRB 990123 with the redshift z ~ 1.6 [2, 56]. 
To explain such a high observable energy release, it is necessary to have a 
strong coUimation, which, in turn, is seriously restricted. Here, we discuss 
various GRB observational data, analyze difficulties and problems of their 
interpretation in the cosmological model, and consider physical restrictions 
on the GRB models. In conclusion, we discuss the problems associated with 
the interpretation of soft gamma repeaters (SGRs) as magnetars and the 
possible relation between short GRBs and SGRs. 

2 Possible energy sources of cosmic gamma- 
ray bursts 

For explaining the GRB radiation, the models of fireball, cannon ball, and 
precessing jet are used. In these models, the key problem associated with 
searching for the possibility of an enormous (10^^10^"^ erg) energy release over 
such a short (0.1100 s) time is avoided. From the various proposed models, 
we select the following ones. 

(i) The coalescence of two neutron stars or a neutron star and a black 
hole of a stellar mass. 

This mechanism was numerically investigated in [86, 87]. The 7 rays are 
produced here due to the {u, v) annihilation, and the energy yield proves 
to be insufficient for explaining the most powerful GRBs even assuming a 
strong coUimation of the GRB radiation. The energy emitted only in the 
GRB 990123 optical afterglow [2, 56] exceeds the total radiation energy yield 
in this model approximately by an order of magnitude. 

(ii) Magnetorotational explosion. 

The magnetorotational explosion proposed in [76] as the cosmological 
GRB model was previously introduced in [15] for explaining supernova ex- 
plosions associated with a stellar core collapse. The 2D numerical calcu- 



3 




Figure 1: Velocity field (on the left) and the specific angular momentum 
v^r (on the right) at the instant of time t = 0.191 s after switching on the 
magnetic field. High-angular-momentum regions are indicated in the right- 
hand figure as the darker ones [71]. 

lations of the explosion of a rotating magnetized gas cloud [5, 6] and the 
calculations of a magnetorotational supernova [71, 7, 8] showed that the effi- 
ciency of transformation of the rotation energy into explosive energy in both 
cases amounts to about 10%. The released energy is sufficient for explaining 
the explosion of a collapsing supernova but proves to be insufficient for cos- 
mological GRB. The results of numerical simulation of a magnetorotational 
explosion are shown in Fig. 1 [71]. 

(iii) Hypernova. 

This model proposed in [76] assumes the possibility of explosion of a 
very powerful supernova. Nowadays, this model is very popular because it 
assumes that traces of supernova outbursts are found in optical afterglows of 
several GRBs [98, 32, 100]. The radiation of a strongly magnetized quickly 
rotating newly born neutron star for the GRB origin was considered in [105]. 
Another model of hypernova assumes the collapse of a massive nucleus, the 
formation of a black hole with a mass M^h ~ 20 Mq surrounded by a massive 
disk, the fast accretion of which results in the GRB phenomenon [60]. This 
model seems to be the most promising now. To the collapse-onset moment, 
massive stars collapsing with the formation of black holes and GRBs were 
most likely [37, 104] in the WolfRayet statea very bright massive compact 
star which lost its hydrogen shell because of the outflow of matter during the 
prior evolution. The SN Ic type of supernova, which is presumably observed 
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at site of GRBs and represents a product of the explosion of a massive star 
devoid of a hydrogen shell, also points to this fact. 

(iv) Magnetized disk around a (Kerr) rotating black hole (RBH). 

This model is based on extracting the RBH rotation energy for the GRB 
production due to its magnetic coupling with the RBH and the surrounding 
accretion disk or torus [106]. 

(v) Dyadosphere. 

In the model [88], the GRB originates from an explosive formation of elec- 
tronpositron pairs and the electromagnetic pulse from an electrically charged 
black hole surrounded by a baryonic remnant. This model is based on the vac- 
uum explosion in the dyadosphere, i.e., the region in which electric- field in- 
tensity supercritical with respect to the formation of e^e pairs is present. The 
key problem here consists in the possibility of formation of such a strongly 
charged BH (sec also [74]). 

(vi) A shock wave beyond the neutron star formed after the supernova 
explosion in a double system was considered as the GRB source in [51]. 

(vii) An exotic model of production of GRBs from superconducting strings 
was proposed in [13]. 

(viii) The phase transition with the formation of a quark (strange) star. 

The energy release at the transition of a neutron star into a stable equi- 
librium state of a quark (strange) star was considered in [14] as the solution 
to the problem of cosmo logical GRB. This model could explain the relation 
of GRB to supernova explosions with the formation of a neutron star and 
the subsequent enormous energy release during its transition into the state 
of a quark star leading to the GRB. An attractive property of this model is 
the possibility of obtaining an arbitrary delay time between the SN explosion 
and the GRB, which is associated with the transition of hadrons in a quark 
state and very strongly depends on various parameters. In this model, the 
delay time can be arbitrarily large so that the majority of SN explosions 
should result in no GRBs according to observations if this time exceeds the 
Hubble lifetime of the Universe. Free quarks should not virtually interact, 
and the density in the case that noninteracting quarks enter an energetically 
favorable state is very uncertain [14]. In Figs. 2 and 3, the models of neutron 
and quark (strange) stars are shown on the massradius plane, and a broad 
spread in theoretical predictions is seen. A greater spread in the properties 
of quark stars in Fig. 3 is obtained in [4] within the framework of the same 
model of quark matter as in Fig. 2, but for a wider variation of parameters. 

In the presence of considerable theoretical problems, only observations 
can settle the question about the possible existence of quark stars. Even in 
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Figure 2: Massradius relation for pure strange quark stars (on the left) and 
hybrid stars (on the right). The G0G4 models for hybrid stars correspond 
to different parameters of the model. H is the pure hadronic star, QC is the 
quarkcore star, and MC is the mixed-core star [101]. 

the presence of theoretical uncertainties, it is impossible to obtain a neu- 
tron (purely hadronic) star with a radius essentially smaller than 10 km. 
The observational discovery of compact stars with such small radii would be 
evidence against their purely hadronic composition. 

3 Statistical properties of cosmic gamma-ray 
bursts 

The conclusion about the cosmological origin of GRBs is based on the anal- 
ysis of their statistical properties and optical-afterglow spectra with strongly 
redshifted lines. 

Statistical arguments for the cosmological origin of GRBs are associated 
with an observed isotropy of the GRB distribution over the sky together 
with a strong deviation of the logiV log 5* distribution (or an equivalent) 
from the Euclidean homogeneous distribution with a slope of 3/2. The ho- 
mogeneous GRB distribution over the celestial sphere was obtained first in 
the experiment KONUS [62] and then confirmed in the experiment BATSE 
[68] in which more than 3000 GRBs were detected in total during nine years 
(19912000) of observation in space. The results obtained in the experiment 
KONUS are shown in Fig. 4; in this experiment, 143 GRBs [63] were detected 
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Figure 3: Massradius dependence for quark stars with A//i = 1.6 and 1. 
Dashed hues represent the models with a weak coupling, which takes place 
for the same choice of renormalization scales [4] . 
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Figure 4: GRBs detected in the experiment KONUS. On the left: the GRB 
position at the celestial sphere represented in the galactic coordinates l^^ and 
b^^ . On the right: log A?" — log 5 distribution [62]. 

in 384 days. 

The authors of [62] assumed that the properties observed for the log N log S 
curve are associated with various selection effects and that the true distri- 
bution of GRB sources in space is homogeneous. The selection effects in 
the experiment KONUS taken into account in [42] led to the average value 
< V/Vmax >= 0.45 ± 0.03; here, 0.5 corresponds to a homogeneous dis- 
tribution of sources. The KONUS data were obtained in 19781980 under 
conditions of a constant background during the flight to Venus. The similar 
analysis of the BATSE data [91] obtained at a circumterrestrial orbit under 
conditions of a variable background gave < V/Vmax >= 0.334 ± 0.008. Fig- 
ures 5 and 6 display the statistical characteristics of the BATSE data [68, 35] 
and the number of GRBs here is more than an order of magnitude higher. 
These two results seem to contradict each other. In the KONUS experiment, 
the sensitivity is approximately 3 times lower than in the BATSE experiment, 
and the deviations from a homogeneous distribution with < V/Vmax >= 0.5 
began much earlier [35]. We show the BATSEPVO joint results in Fig. 7 
from [35] . The PVO data on the brightest GRBs show good spatial unifor- 
mity for GRBs with a slope of 3/2. In [92], M. Schmidt made a detailed 
statistical analysis and processed the BATSE data by dividing the GRBs 
into four classes according to the degree of spectrum hardness; he calculated 
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Table 1: 1. Dependence of < V/Vmax > on the degree of hardness for 1391 
GRB [92] 



number 


< Q;23 >obs 


< V/Vrnax >obs 


< Q;23 >corr 


^ ^/^max corr 


348 


-2.55 


0.468 ±0.017 


-2.33 


0.421 


348 


-1.84 


0.309 ±0.016 


-1.79 


0.325 


347 


-1.47 


0.299 ±0.016 


-1.47 


0.344 


318 


-1.01 


0.270 ± 0.015 


-1.10 


0.256 



< V/Vmax > for each class individually. The results of this investigation are 
listed in Table 1, where determines the spectrum slope obtained from 
the counts in the second and third BATSE channels with the energy ranges 
of 100300 keV and 50100 keV, respectively. In the column obs, we hst the 
observables corrected taking into account statistical errors in count peaks 
and, in the column corr, we give the same data assuming the existence of a 
luminositydegree-of-hardness correlation. 

In the cosmological model, it was possible to expect a lower value of 

< V/Vmax > for softer GRBs in the case of a homogeneous set because, with 
an increase in the redshift, the spectrum becomes softer. The result proved 
to be the oppositesoft sources had a greater value of < V/Vmax > than the 
hard ones: 0.47 and 0.27, respectively. 

In [92], it was assumed that the increase in the GRB power with the degree 
of spectrum hardness is so significant that it outweighs the opposite effect 
caused by the redshift in a homogeneous sample of objects. Another expla- 
nation can be that the sample of softer GRBs is more complete because of 
selection effects. Figures 8 and 9 illustrate a possible key role of the selection 
effects (incompleteness of the data set and statistical errors in estimates of 
luminosity in the presence of threshold). The incompleteness of data affects 
the distribution of such well-investigated objects as G stars of the Sun type. 
It is possible to expect still greater effects when detecting such short tran- 
sients as GRBs. The comparison of average values of V/V^ax as functions of 
the detection threshold for G stars and GRBs in Fig. 8 from [41] shows their 
qualitative agreement. 

The effect of statistical errors in the presence of a detection threshold is 

investigated in [17]. The V/Vmax distributions with average statistical errors 
at the level of ten thresholds are shown in Fig. 9a for the normal distribution 
of counts and in Fig. 9b for the similar distribution over their logarithms. 
The distribution in Fig. 9a is similar to the BATSE distribution in Fig. 6. 
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Figure 5: Distribution of 1121 GRBs over the celestial sphere in galactic coor- 
dinates detected by BATSE over three years. No clusterization or anisotropy 
is observed [35]. 

4 Optical afterglows and redshifts 

The x-ray afterglow discovered by the Italian Beppo-SAX satellite made it 
possible to identify GRBs optically and to obtain their spectra. These spectra 
showed the presence of a large redshift z, reaching 4.5, which points to the 
cosmological origin of GRBs and an enormous energy release. In most cases, 
the redshifts were measured for very weak parent galaxies. The list of certain 
redshifts is given in Table 2, where the redshift data taken from [34] are 
supplemented with the total radiated GRB energy [19]. In the table, we 
listed the numbers of triggers and total energy fluxes of GRBs from the 4B 
catalogue [75] and total energy fluxes from other sources. The location of 
optical afterglows in parent galaxies is given in [26]. The enormous energy 
yield in a short time (from 0.1 s to several hundreds of seconds) is the primary 
difficulty in the cosmological interpretation. In certain cases, for example, 
in GRB 990123, the burst isotropic energy (2.3 ■ 10^^ erg) exceeds the rest 
energy of the Sun. All of the mechanisms discussed above, with the exception 
of the Dyadosphere model of Ruffini and collaborators, could give an energy 
release less than one percent of this quantity. In the case of Dyadosophere 
the Blackholic energy could fulfill the above requirement 
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Figure 6: (a) Integrated distribution of the GRB number as a function of 
the maximal counting rate for 140 GRBs. For a homogeneous distribution of 
sources, the power law with an exponent 3/2 is expected. The total detection 
rate for bursts amounts to ~ 800 events per year, (b) The distribution of 
V/Vmax is given for 140 GRBs. An average V/Vmax is equal to 0.348±0.024 
[68]. 
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Figure 7: Distribution logA^ — logP for a combined set of BATSEPVO 
data. The distributions match well in the dataoverlap region. It can be seen 
that the PVO data containing a larger number of powerful bursts than the 
BATSE data and obtained for longer time well follow the power law with an 
exponent 3/2 for powerful bursts [35]. 
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Figure 8: Average value of V/Vmax as a function of the minimal brightness 
of the objects taken into account for G stars (closed circles) and GRBs from 
BATSE 2B catalogue (crosses with the indication of the la-error interval). 
For G stars, the error intervals are less than the circle sizes. The crosses in 
circles give a value of average V/Vmax from the BATSE IB catalogue. The 
hmit of the GRB-detection estabhshed with BATSE corresponding to 5.5(7. 
For G stars, the hmiting value m„ = 10 corresponding to S/N = 5.5 was used 
in this analyses. Note that < V/Vmax >■ is close to the uniform value 0.50 
both for bright G stars and for bright GRBs. By taking weaker sources into 
account, < V/Vmax > deviates in the lower-value side. For G stars, this fact 
is associated with an incompleteness of the catalogue of weak sources [41]. 
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Figure 9: (a) Curve [log N — log C (max) / C (thr)] in the presence of statistical 
errors normally distributed with an average error Ai in detection-threshold 
units; 1, is the straight line with a slope of 3/2 corresponding to Ai = 0; 2, 
is the curve with Ai = 1; and 3, is the curve with Ai = 10. (b) The same 
as in Fig. 9a for the normal logarithmic distribution of errors; A determines 
the average error as a spread of the threshold- number logarithm; 1, is the 
straight line with a slope of 3/2 corresponding to A = 0; 2, is the curve 
with A = 1; C(max) is the peak burst intensity; and C(thr) is the value of 
detection threshold [17]. 
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Table 2: Redshifts of parent galaxies (Z) and GRB total energy fluencies (F), 
June 2001 



Trigger 


GRB 


R mag 


Z 


Type-^ 




Ref. 


number 










erg/cm^ 
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^ o n 
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d7U7 


nor\ /I o c 

980425 


14 


n nno c 

(J.UU85 


a,e 


A A 1 n — 6 

4.4 • ID 


[36] 


d7d4 


98U519 


25.2 
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9oUd1o 


24. U 
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e 




h 1 /II 
[114j 


0591 


9»L)7L)3 


22. D 


n n/^c 

0.966 
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1^ A O 

24. o 
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'70 A O 

7343 


99U123 


oo n 

23.9 


1 £?nn 

1.600 


a,e 


CI 1 n — 4 

5.1 • 10 


[52] 


7457 


99U3L)o 


>2».5 






in 1 n — 5 /o \ A\ 

1.9 -10 (3 + 4j 




7549 
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1 d n 
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99U7L)5 


22.0 


n or' 

0.86 
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1 n — 5 

~ 3 • 10 


[61] 




990712 


21.8 


0.4331 
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QQ1 90^ 




U. 1 uoo 
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in-4 

~ lU 


1331 
[44j 


7906 


991216 


24.85 


1.02 
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2.1 • 10-4(3+4) 




7975 


000131 


>25.7 


4.50 


b 


~ 10-5 


m 




000214 




0.37-0.47 


X 


~ 2 • 10-5 


m 




000301C 


28.0 


2.0335 


a 


~ 4 • 10-6 


m 




000418 


23.9 


1.1185 


e 


1.3- 10-5 


m 




000630 


26.7 






2 • 10-6 


m 




000911 


25.0 


1.0585 


e 


5 • 10-6 


m 




000926 


23.9 


2.0369 


a 


2.2 • 10-5 


m 




010222 


>24 


1.477 


a 


The brightest on 
BeppoSAX 


m 



° e is the emission line, a is the absorption line, b is the continuum break, and x 

are the x-rays. 

Relation of this galaxy with SN or GRB is not quite reliable. 
Relation of optical transient with this GRB is unreliable. 



^ Number of B ATSE channel with a maximal flux is indicated in parentheses [75] ; 
otherwise, the estimates of the total energy flux are given from other references. 
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4.1 Collimation 



To avoid an enormous production of energy, it is assumed that the GRB 
radiation is colhmated. In the cannon-ball model [32], it is assumed that a 
whole radiating object moves relativistically with a large factor F Ri 10^ — 10^ 
that leads to a coUimation solid angle Q 10~^ — 10~^. The analysis of the 
GRB collimation was made in [85]. A strong collimation-angle restriction 
follows from the analysis of the probability of the appearance of orphan 
optical afterglows, which are most likely poorly collimated or not collimated 
at all. The absence of certain rapidly alternating orphan optical objects 
results in the following restrictions. In the case of isotropic GRBs, it was 
expected to detect ~ 0.2 orphan afterglows so that the absence of such objects 
assumes that Qopt/^-y 100, which suffices for ruling out an especially strong 
collimation. The similar investigation using the data on afterglows in the 
radio-frequency band and the variable radio sources made in [79] leads to 
the restriction on the coUimation angle 9j > 5°. Because the afterglows in 
the radio-frequency band occur at the nonrelativistic stage of evolution of a 
GRB remnant, it is expected to be isotropic, and the restrictions following 
from the absence of radio orphans at Vtr/Vt^ establish a limit on the radiation 
collimation Q-y of the GRB itself. 

The comparison of redshifts and total energy fluxes from Table 2 shows 
no correlation between the distance and observable flux (sec Fig. 10). This 
fact is usually explained by a strong collimation leading to a wide spread in 
observed fluxes due to observations of a narrow beam under different angles 
of view. If the collimation were associated with the relativistic motion of a 
source [32] , the strong correlation between the observed radiated energy and 
the GRB duration would be detected in this case: powerful GRBs should 
have been shorter. The absence of such a correlation allows one to rule out 
the models based on the ultrarelativistic motion of sources. 

4.2 Prompt Optical Afterglows 

The optical afterglow of GRB 990123 was successfully observed 22 s after 
the beginning of its detection in the gamma region [3, 2]. GRB 990123 was 
detected by BATSE in 1999 on January 23.407594. The burst was strong with 
duration >100 s, involving several peaks (see Fig. 11), and revealed a strong 
spectrum evolution. The T50 and T90 durations at the 50 and 90% energy- 
flux levels were equal to 29.82 (±0.10) and 63.30 (±0.26) s, respectively. The 
peak of optical brightness at a level of 8.95™ was achieved within 30 s after 
the burst onset, and the brightness decreased down to 14.5"* within 95 s. 
Thus, the optical-radiation peak almost coincided, with a certain delay, with 
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Figure 10: Total observed energy fluence F as a function of the redshift z for 
GRBs from Table 2 [20]. 
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the gamma ray peak. The observed optical luminosity, which was related 
to the redshift z — 1.61, achieved Lop* fa 4 • 10^^ erg/s, which is brighter by 
almost five orders of magnitude than the optical luminosity of any supernova. 
The energy release in the direct optical radiation reaches 10^^ erg, and an 
isotropic gamma-ray flux exceeds 2.3 ■ 10^^ erg, which is more than the Sun 
rest energy [2, 56]. Another bright afterglow was observed for GRB 021004 
(15"^, z = 2.3), GRB 030329 (12.4"*, z = 0.168), and GRB 030418 (16.9"*). 
In brackets, we list the brightest visual magnitudes and redshifts. 

The GRB 030329 afterglow observed at many observatories is most in- 
teresting (see, for example, [89]). Features in the afterglow spectrum showed 
the presence of a supernova [100]. Figure 12 shows the .-ray observations 
for GRB 030329 in the KONUS WIND experiment, and Fig. 13 taken from 
[82] shows the hght curve for the optical afterglow obtained at the Crimean 
observatory. 

5 GRB and supernovae 

The observational indications of the relation to supernovae are found for the 
following GRBs: GRB 980425 {z = 0.0085, 40 Mpc), GRB 980326 {z^l), 
GRB 011121 {z = 0.365), GRB 020405 {z = 0.695), and GRB 030329 {z = 
0.169). The light curves for certain optical afterglows have a red bump for 
1575 days, which may be induced by an accompanying supernova explosion 
(see [98], Fig. 16 from [27], and Figs. 18 and 19 from [43]). The indications 
of the relation between GRBs and supernova explosions were obtained in the 
observation of GRB 980425, which coincided in the sky with SN 1998bw. 
This supernova proved to be abnormally bright in comparison with other 
SNe (such as SN Ib/c) and had peculiar spectra and unusually high radio- 
luminosity at initial stages [78]. It exploded in a nearby galaxy with a much 
smaller redshift than those observed in the optical afterglows of other GRBs 
from Table 2, which points to a low energy release in GRB 980425. It is 
impossible to rule out that the SN and GRB approximate coincidence in 
time and projection to the celestial sphere is accidental, and the distances to 
them are very different. Nevertheless, even given the obvious peculiarity of 
nearby GRB 980425 with its relation to SN 1998bw, the idea of GRB origin 
in the explosion of the very bright peculiar supernova SN Ic (hypernova) 
became very popular. The features of light curves of optical afterglows and 
red bumps associated with their nonmonotonic behavior were interpreted 
as the accompanying supernovae, which would otherwise remain unnoticed 
because of the very large distances to a host galaxy. 

Another piece of evidence in favor of the GRBSN relation is obtained 
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Figure 11: Thin lines reproduce the gamma-ray profile from the BATSE 
data with the resolution of 1024 ms in various energy channels. The thick 
line connects some first dots of ROTSE optical observations, which began 
within 22 s after the detection. This burst is the brightest one in the optical 
range [82]. 
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Figure 12: GRB 030329 in gamma rays: KONUSWIND (UT)11:37:29. The 
total energy fluence is equal to 1.2 x 10""^ erg/cm^, the duration is 50 s, and 
the maximal flux is 2.5 x 10~^ erg/(cm^s); one of the brightest bursts in the 
gamma region [82]. 
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Figure 13: Hour-scale variation of the GRB 030329 optical afterglow [82]. 
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Figure 14: On the left: the GRB 030329 afterglow spectrum. This spectrum 
covering a wide region is obtained within 0.5 days after the GRB. Closed 
circles correspond to the measurements obtained nearby the nominal time, 
and open circles correspond to the extrapolation to the nominal time assum- 
ing the evolution in a constant-density medium. The simple fit of this broad 
spectrum gives the following parameters: the synchrotron self-absorption fre- 
quency Ua ~ 25 GHz, the spectral-peak frequency Um ~ 1270 GHz, the cool- 
ing frequency Uc ~ 6.2 x 10-^^ Hz, the spectral-peak flux ~ 65 mJy, and the 
electron-energy index p ~ 2. The following physical parameters correspond 
to these values: the explosion energy £^ ~ 5.7 x 10^^ erg, the surrounding-gas 
density n ~ 5.5 atom cm^, the electron-energy fraction eg ~ 0.16, and the 
magnetic-field-energy fraction ~ 0.012. On the right: the light curve of 
the GRB 030329 optical afterglow in the R region of the spectrum beginning 
approximately after one day from the GRB moment. The observation errors 
are less than the data point sizes [83] . 
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Figure 15: Spectrum of optical afterglow obtained with the telescope Keck I 
using the spectrometer with an effective resolution of 4.2 A. The exposition in 
the afterglow observation amounted to 600 s. The identified narrow emission 
lines [O 11], HjS, [O 111], and the lines of Mg 11 absorption correspond to 
an average redshift z = 0.169 ± 0.001, which makes GRB 030329 to be a 
cosmological GRB with the smallest redshift. The identified emission hues 
are typical for star-formation regions in galaxies, while the absorption lines 
are caused by gas in the galaxy disk. The identified line Ca II at z ^ is 
hypothetically associated with the absorption in a cloud inside our Galaxy 
[83]. 
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Figure 16: (a) Light curves of afterglow and (b) red bumps observable at 
an intermediate time in GRB 011121. The triangles correspond to the pho- 
tometry with the Hubble space telescope (HST) with the F555W, F702W, 
F814W, and F850LP filters (everything with corrections on an estimate of 
contribution from the host galaxy), and the rhombs represent the earth ob- 
servations. The dashed line describes an assumed behavior of GRB optical 
afterglow, and the dotted line sets an expected flux from a sample supernova 
for the GRB 011121 redshift, which approximately corresponds to the obser- 
vational data with taking into account the 55% absorption. The solid line 
corresponds to the sum of radiations from SN and GRB afterglow [27]. 
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from the detailed photometric and spectroscopic investigation of the bright 
optical afterglow of GRB 030329. At a certain stage, it was noted [100, 43] 
that the spectra of afterglow and peculiar supernova SN 1998bw are similar; 
this allowed the authors to declare the discovery of SN 2003dh associated 
with this GRB. 

We give the comparison between the spectra of SN 2003dh and SN 1998bw 
in Figs. 17, 18, and 19. In Fig. 18, the upper spectrum is very well approxi- 
mated by a power law, which is usually observed in the spectra of afterglows. 
The intermediate spectrum shows an appreciable difference from the power 
law similar to that observed in SN 1998bw at the same phase. At the lower 
spectrum, where SN 2003dh dominates, the supernova features are clearly 
seen. For comparison, we show the spectrum of SN 1998bw at the 33d 
day after the flare (the dashed hue) shifted according to the GRB 030329 
redshift. All the SN 2003dh spectra are presented in observed wavelengths 
without taking into account corrections on redshift. The lines of absorption 
in the atmosphere were also omitted. In the range of wavelengths longer than 
9000 A, the spectrum is strongly polluted with lines of the night airglow, but 
the wide detail around 10 000 Acan be associated with a supernova. At all 
stages, the emission lines [O II]A3727, H^, [O III]A4959 and A5007, and H« 
are identified as the lines most likely belonging the host galaxy. At the last 
stage (May 1), the following lines also were identified: [Ne III]A3869, H5, H7, 
[N II]A6583, and the superimposed lines He IA3889 + H8 [Ne III]A3968 + 
H^. Strong Balmer lines point to the fact that the absorption in the host 
galaxy is weak. The estimate of metallicity based on fiuxes in the lines [O 
II], [O III], and Hjs gives [O/H] = 1.0. The estimate of a star-formation 
rate as 0.2 M© yr— 1 is obtained based on the lines [O II] or Taking into 
account the 3a upper hmit R > 22.5"* obtained for the emission of the host 
galaxy when studying archive data, it was concluded that equivalent widths 
of emission lines are very large. Thus, the host galaxy is the dwarf with a 
low metallicity and an active star formation that qualitatively coincides with 
the host galaxy of GRB 980425 and SN 1998bw. 

In Fig. 19, the solid lines set the SN 2003dh spectra obtained as a result 
of the processing described in [43]. The dashed lines set the SN 1998bw 
spectra observed at similar epochs. We give all the spectra in the restframe 
wavelengths and marked by indicating days from the GRB-fiare moment with 
taking into account the cosmological time extension (1 + z). The relative 
values of spectra along vertical are arbitrary. The SN 2003dh spectrum was 
divided into 15-Aregions, the lines of the host-galaxy radiation were removed, 
and an interpolation was made where it was necessary in the region of strong 
night-sky absorption lines at 6800 A, 72007400 A, and 7600 A. The strong 
absorption "edge" at 6150 A(observed at 7200 A) is likely associated with 
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Figure 17: Spectrum of SN 2003dh obtained on April 8 after subtracting 
of the smoothed spectrum of April 1. This residual spectrum shows wide 
bumps in the regions of approximately SOOOAand 4200 A(in the rest frame), 
which are similar to those observed in the spectrum of a peculiar SN 1998bw 
Ic type within one week before the brightness peak [78]. The coincidence is 
not so good compared with the spectrum of SN 1998bw in the brightness 
peak especially for the red spectral edge [100]. 
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Figure 18: Evolution of the spectrum of total optical radiation for the GRB 
030329 afterglow, the SN 2003dh supernova related to it, and their parent 
galaxy [43]. 
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Figure 19: Comparison between the evolution of SN 2003dh and SN 1998bw 
spectra [43]. 



28 



the night sky absorption. When processing the spectra, it was assumed that 
these are only the host-galaxy emission lines that were visible and subtracted. 
The spectrum was fit in the form of the sum of a power part {fx ~ A_(^+2)) 
and the scaled spectrum of the sample SN 1998bw. Using the least square 
technique, three parameters have been obtained: the spectral index f3, the 
afterglow amplitude, and the supernova amplitude. In most cases, the best 
index was obtained equal to /3 ~ —1.2 ± 0.05, which was used everywhere; 
however, this index and the sample spectrum shape only slightly affected the 
estimate of contribution of supernova in the spectrum shape. The spectral- 
peak wavelength redshifts in time for both supernova, and the average rate 
of such a shift amounts to ~ 25Aper day for SN 2003dh that is similar to the 
early evolution of the SN 1998bw spectrum. In [43], an increase in opacity 
on absorption is assumed as the cause of such a shift for the motion in bluer 
region from 4900 A(rest-frame wavelength). 

Figure 20 shows the earliest optical-afterglow spectra of GRB 030329 ob- 
served in [99]. From the comparison between Figs. 18 and 20, it follows that, 
at early moments, for example, within 1012 h after the GRB, the spectral 
shape was far from that described by the power law to which it evolved within 
several days to April 3. In Fig. 20, it could be interpreted as the spectrum 
of an initial strong outburst of SN 2003dh similar to a short initial outburst 
observed on the light curve for SN 1993 J shown in Fig. 21. The model of 
such an outburst was calculated in [95]. 

Despite such striking coincidences, the relation between GRB and SN still 
cannot be considered as firmly established. There are spectral distinctions, 
time scales differ from standard values, and the spectra of many GRB after- 
glows involve no indications of the presence of supernovae. In addition, in 
a unique case when the coinciding GRB and SN were discovered indepen- 
dently, a quite ordinary observable GRB 980425 appeared together with SN 
1998bw proved to be very abnormal, being three orders of magnitude weaker 
than other cosmological GRBs with known redshifts of lines in afterglows 
(see Table 2). 

6 Polarization of gamma rays in GRB 

The authors of [30] reported on the discovery of the linear polarization of . 
rays from GRB 021206. The authors interpreted this fact as a synchrotron 
radiation of relativistic electrons in a strong magnetic field. Measurements 
were carried out at the Rcuvcn Ramaty High Energy Solar Spectroscopic Im- 
ager (RHESSI). The RHESSI has a set of nine coaxial germanium detectors of 
a large fiducial volume (300 cm3) with a high spectral resolution for studying 
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Figure 20: Smoothed spectra of GRB 030329 optical afterglow in rest-frame 
wavelengths obtained with the spectral resolution of 12 A by the 6-m tele- 
scope SAO. The spectra obtained at the time moments of 11.5, 11.3, and 
10.8 h after the GRB are shifted along the fx vertical scale by +0.2 • 10~^^, 
+0.6 • 10-^5 and +1.2 • 10"^^ respectively, [99]. 
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Figure 21: Light curves calculated using the model with a shell mass 
O.89M0. Bolometric luminosity measurements for SN 1993J are given 
closed circles [95]. 
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X and 7 rays from the Sun in the energy range of 3 keV17 MeV. The RHESSI 
has a high angular resolution (2") with a field of view of ~1° of its optics; 
however, flat focal detectors are unscreened and cover the entire sky. Because 
of this fact, the RHESSI very frequently detects GRBs in the focal plane of 
detectors, although the probability to fall into the image-construction field is 
low. These observations allow the authors to obtain highly resolved spectra, 
while detection times and energies of individual photons are potentially suit- 
able for measuring the polarization. The RHESSI is not optimized for use as 
a 7-polarimetcr, but some features of its design make it the most sensitive 
of modern devices for astrophysical measurements of the 7-ray polarization 
[30]. 

For the soft 7 rays with ~0. 15-2.0 MeV, the basic mechanism of interac- 
tion of photons with RHESSI detectors is the Compton scattering. A small 
fraction of incident photons undergoes the single scattering in one of detectors 
before being scattered or photoabsorbed in another detector. The probabil- 
ity of this sequence of events depends on the polarization of the incident 
7 ray. The linearly polarized 7-rays scatter predominantly in the direction 
perpendicular to their polarization vector. In the RHESSI, this property of 
scattering can be used for measuring the polarization of astrophysical sources. 
The device sensitivity to the polarization depends on an effective scattering 
area and an average value of polarimetric modulation depth ij,{9,E), which 
is equal to a maximal spread in the probability of azimuthal scattering of 
polarized photons. This depth is equal to /i — {da± — da\\)/{da± + da\\), 
where da^, and da\\ are the KleinNishina differential cross sections for the 
Compton scattering in perpendicular and parallel directions relative to the 
polarization direction, respectively. It is a function of the incident-photon 
energy E"^ and the Compton scattering angle 9 between the direction of inci- 
dent and scattered photons. For a source with the count number S and the 
degree of polarization H^, the distribution of expected azimuthal scattering 
angles is equal to dS/d(j) = (S'/27r)[l — fimXis cos{2 {(f) — r]))] , where (f) is the 
azimuthal scattering angle, rj is the polarization- vector direction, and Hm is 
the average value of polarimetric modulation depth of the device. Although 
the RHESSI has a small effective area (~20 cm^) for the events in which the 
scattering between detectors takes place, it has a very large modulation in 
the region of 0.152.0 McV obtained in the Monte Carlo calculations. Ac- 
cording to [30], the results of processing of observational data shown in Fig. 
22 convincingly point to the discovery of a strong polarization for the GRB 
direct 7 rays. The authors of [30] assumed that this result specifies the im- 
portant role of the magnetic field in the GRB-explosion mechanism. It agrees 
with the GRB model in the form of a fireball inflated by the magnetic field. 
This field could be generated due to the rotational energy of an accretion 
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disk differentially rotating around a central compact object, or the rotation 
energy of the Kerr black hole pierced by magnetic field hnes, or the rota- 
tional energy of a strongly magnetized neutron star [105]. Another possible 
model is the generation of a large-scale magnetic field in a flux beyond the 
shock- wave front. 

Data [30] were analyzed again in [90], where it was concluded that they 
give no hint of the detection of a polarized signal (see Fig. 23). First, it 
was found that the number of polarization events used for the measurement, 
where the photon passed from one detector to another as a result of scattering 
and was detected in both, was approximately ten times lower than that 
estimated in [30] 830 ± 150 instead of 9840 ± 96). As a result, the signal- 
to-noise ratio proves to be too low for detecting a possible polarization even 
at the level of 100%. In [90], the authors used another method of data 
processing for detecting the polarization effect and revealed no indications 
of the polarization in the GRB 021206 radiation. It was shown that the 
detected signal can correspond to an unpolarized source, and the polarization 
observed in [30] is associated with neglecting a systematic error in the zero 
light curve used in processing. Because of the low signal-to-noise ratio in the 
RHESSI data, the analysis made in [90] by a method which can be employed 
only for the analysis of Poisson noise, demonstrated that the data could be 
consistent with any polarization up to 100%. This points to the fact that it is 
impossible to obtain estimates of polarization of the GRB 021206 radiation 
from the RHESSI observations. 

In reply to [90], the authors published a new paper rejecting this crit- 
icism. It was noted that [90] "raises many important problems associated 
with processing the observations, which will be considered in a separate pa- 
per being prepared for publication. However, it should be noted that the 
limit for the detected degree of polarization obtained in [90] with the use 
of the new method of processing developed there is strongly overestimated. 
Although we would like to learn in more detail the new statistical technique 
proposed in [90] (especially with the substantiation of the assumption that 
the obtained degree of polarization is independent of the instrumental re- 
sponse), we come to an inevitable conclusion that there is a serious gap in 
their statistical method. In addition, the authors of [90] declared that their 
analysis is insensitive to an arbitrary level of polarization. Hence, it does 
not contradict that level of polarization, which was presented in our orig- 
inal paper." The investigation performed in [113] confirmed the conclusion 
of [90] that it is impossible to determine the degree of polarization of GRB 
radiation from the available data. It is possible that the further processing 
method development and new observations are necessary for overcoming this 
contradiction; however, the majority is on the side of skeptics. 
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Figure 22: On the left: RHESSI light curves (total number of counts) in 
three energy intervals for GRB 21206. This GRB is localized in the sky at 
18° from the Sun, ruling out the possible search for the optical afterglow; 
however, its brightness, duration, and proximity to the RHESSI rotation 
axis make it an ideal candidate for searching the polarization. The shaded 
region shows the signal with a 5-s averaging time for the analysis of polariza- 
tion. On the right: the azimuthal distribution of scatterings for the RHESSI 
data corrected for the satellite rotation. The angular interval of 0°-180° is 
divided into subintervals of 15° in which the counts represented here were 
gathered twice for the sake of clarity. The upper graph shows the unpro- 
cessed measured distribution (crosses) together with the result of the Monte 
Carlo simulation for a simulated distribution of an unpolarized source with 
a constant flux (rhombs). The lower graph shows the RHESSI data after 
subtracting the simulated distribution. The obtained residual distribution is 
incompatible with the unpolarized source (the dashed line) at the confidence 
level higher than 5.7a. The solid line shows the result of the best flt of the 
modulation curve corresponding to the degree of polarization of (80 ± 20)% 
[30]. 
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Figure 23: Here, ^ = ^„=0 correspond to the direction towards the celes- 
tial north with an increase in the angle when rotating the detector base from 
the north to the east (clockwise facing the Sun). The GRB 021206 position 
in all the figures is shown by the vertical dotted line, (a) The counting rate 
S of soHtary events as a function of 9, where 9 is the angle between the ce- 
lestial north and the center of the detector detecting a photon in the interval 
of 0-360°. S{6) is incompatible with a constant value, C,l is indicated; the 
observation probability of such a variability is lower than 10^^'^'^. (b) The 
counting ratei^ of double events as a function of 6; here, both counts are 
indicated for each event located in the corresponding intervals of 9. The ob- 
served variability is incompatible with the constancy of D[9), is indicated; 
the observation probability of such a variability is lower than 10~^°. (c) The 
ratio S/D is compatible with constant with the probability equal to 0.54. 

(d) The counting rate S{9n) of solitary events for the polarization analysis 
in the interval of 0-180°. The counting rate is incompatible with constant. 

(e) The counting rate D{9n) of double events for the polarization analysis in 
the interval of 0-180° is well compatible with the corresponding graph from 
[30] without subtracting the Monte Carlo zero light curve. An essential vari- 
ability on 9n is seen, (f) The ratio R = D{9n) / S{6n)- The polarized photons 
would scatter in a sky predominate direction with a higher ratio of double- 
to-solitary counts than in the perpendicular direction leading to a variability 
in R{9n). From the observations, it follows that R{9n) is compatible with 
constant; i.e., no indications of polarization were found [90]. 
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Nevertheless, even without the assurance in rehabihty of measurements of 
the polarization for the direct GRB radiation, certain theoretical models were 
developed for explaining this phenomenon with the use of both synchrotron 
radiation [38] and inverse Compton scattering [57]. In [38], the degree of 
polarization P ~ 43 — 61% was found for an ordered transverse magnetic 
field -Bord- At the same time, the field Bi^ generated on the shock-wave front, 
which is chaotic but located completely in the shock-wave plane, can provide 
the degree of polarization up to P < 38 — 54% in a solitary pulse on the GRB 
light curve. However, it is expected that the polarization in the integrated 
radiation of several pulses measured for GRB 021206 decreases approximately 
by half. The magnetic field Py normal to the shock-wave front can provide 
the degree of polarization P ~ 35—62% for the total radiation of many pulses. 
However, the measurements of the polarization in GRB afterglows assume a 
more isotropic configuration of the field generated by the shock wave that 
should reduce P ~ 23 times. Therefore, the ordered field Pord arising in a 
source can explain the observed polarization most naturally, while Py is less 
probable, and P^ is least expected. 

In [57], the inverse Compton scattered radiation of a fireball was con- 
sidered as the directed jet with an opening angle comparable to or larger 
than the collimation angle associated with the relativistic contraction of the 
beam. The degree of linear polarization as a function of the jet-opening angle 
and the direction to observer were numerically calculated for the Compton 
upward scattering along the flux for a large-scale motion (see Fig. 24). Note 
that the gamma factor of the large-scale motion cannot be large due to the 
absence of apparent inverse dependence of the GRB duration on its power, 
which is inherent to the model with a large-scale relativistic motion (see 
section 4.1). 

7 Observation of x-ray echo 

The observations of GRB 031203 at the Newton x-ray observatory (XMM- 
Newton) began on December 4, 2003 at 04:09:29 UT and lasted 58 211 s. 
This GRB was detected for the first time by the IBIS device on the satellite 
INTEGRAL 2003-12-03 at 22:01:28 UT. The analysis of the first observations 
by the Newton telescope found a diffusive x-ray halo around the afterglow 
site. All three chambers of the EPIC device observed this halo, and it is not 
associated with the scattering of optical or x-ray photons inside the device. 
The halo shape was almost a perfect circle, the radius of which increased in 
time. This recalls the behavior of the light echo going from x-ray photons 
scattered by the dust of a cloud located at a distance of about 700 pc from 
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Figure 24: On the left: the degree of polarization as a function of the obser- 
vation angle 6obs in units of 1 /F for a homogeneous jet with sharp boundaries. 
Various lines represent the polarization for jets with various opening angles. 
A bolder line designates the region where the scattering efficiency is higher 
than 2.5%. On the right: the same picture for a jet with the Gaussian dis- 
tribution. Because the Lorentz factor is variable in this jet, the observation 
angle 9obs is indicated along the axis x in units of I/Fq, i.e., in units of the 
Lorentz factor along the jet axis [57]. 
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the observer [59]. GRB 031203 has galactic coordinates 1 = 255.74° and 
b = 4.80° in the direction of which the Gum Nebula, other nebulas, and 
infrared sources are located. The distance to the scattering cloud points to 
its location inside our Galaxy. The x-ray spectrum of GRB 31203 is well 
approximated by the power law with a photon index of ~ 1.7. As expected, 
the scattered x-rays have a softer spectrum with a photon index of ~ 3 [59]. 
The subsequent analysis shows [107] that the halo appears in the form of 
concentric circular rings with the center at the GRB localization (see Fig. 
25). The radii of these rings grow in time as ~ t^/^, which is compatible with 
a small angle of scattering of x-ray photons on a large dust column along the 
direction of observation of the cosmological GRB (Fig. 26). The appearance 
of two rings is associated with two different dust layers in the Galaxy located 
at distances of 880 and 1390 pc, which agrees with known features of the 
Galaxy structure. The halo brightness points to an initial momentum of the 
soft X rays emitted simultaneously with the observed GRB. 

We note that the same ~ t^/^ law for an increase of the ring radius in 
time takes place for an arbitrary variant when the ring radius {x) is much 
less than the distance from a dust cloud to both the observer (di) and the 
GRB source (^2)- Always assuming that x <^ di, we obtain 



X 



= diV2 



- + ^ + 1 - 

di di \ 




1/2 

(8) 



Here, ^2 is arbitrary, and the time is counted from the moment of detection 
of the direct GRB x-ray (F-ray) pulse. When the source is inside the cloud 
((^2 = 0), we have from Eq. (8) that 

X ^ ct. (9) 

In the case oi d2 x, we have 




i.e., the same square-root time dependence. Here, the only requirement is 
large distances di and ^2 in comparison with x, which can be fulfilled even 
for the GRBs located at noncosmological distances. Thus, an increase in the 
radius x{t) according to the square- root law is compatible with the cosmo- 
logical GRB nature, but cannot be considered as proof of this nature. 
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Figure 25: Series of images of the region in ten angular minutes around GRB 
031203 in the energy region of 0.72.5 keV. The observational data are divided 
into ten time intervals of 5780 s. The first six images are smoothed with the 
use of a Gaussian nucleus of six angular seconds in size. The three brightest 
point sources (SI, S2, and S3) are marked [107]. 
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Figure 26: Time increase in the radius of two rings around GRB 031203. 
The radius of each ring was measured from the local peak along the radius. 
In both cases, the distribution is well approximated by the function 9 oc 
{t — toY^"^. A segment length used in measurements of profile along radius 
was taken as an error [107]. 
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8 Reaction of a dense molecular cloud to a 
nearby cosmological GRB 



8.1 Optical afterglow from interstellar-gas reradiation 

As was noted for the first time in [76] (see also [81]), the GRB-afterglow 
properties are better explained under the assumption that GRB sources are 
located in the star-formation regions with a high gas and dust density. The 
interaction of a powerful GRB pulse with the surrounding gas of a density 
n = 10^ — 10^ ~^ results in the occurrence of a special-form optical afterglow 
the duration of which can reach ten years. The light curve and the spectrum 
of such an afterglow were calculated in [25] (see Fig. 27). It was shown 
that the duration and spectrum of the GRB optical afterglow depend on 
the GRB location with respect to the dense cloud because of the interaction 
of 7 rays with the dense interstellar medium. In the case of an anisotropic 
GRB, they depend also on the orientation of the observer and the cloud with 
respect to the GRB axis. During the explosion in a homogeneous cloud, the 
optical afterglow can endure for ten years. The GRB optical afterglow can be 
distinguished from a supernova radiation with the same energy yield by the 
light-curve and spectral features. The Ha and Hf^ emission lines are strongest 
in the optical afterglow. The discovery of the similar afterglow would give 
important information on the GRB-source properties and the surrounding 
medium. 

More detailed 2D calculations of the GRB afterglow after the explosion 
in the star-formation region nearby a cold molecular cloud were made in 
[9, 10]. Optical light curves qualitatively coincide with those in Fig. 27; 
however, the large variety of spectra and light curves takes place depending on 
initial distributions of density in cloud, the anisotropy of density distribution, 
and GRB radiation. For solving the equations of hydrodynamics with the 
radiation heating and cooling in the complete-transparency approximation, 
the PPM method was used. As in [25], the GRB propagation through a cloud 
was approximated by a very narrow spherical or directed gamma-photon wave 
leaving behind a completely ionized gas heated by the Compton interaction. 
When calculating the heating, the KleinNishina corrections to the Thomson 
cross section [11] were taken into account. The temperature distribution 
after an isotropic GRB explosion in the center of a homogeneous spherical 
cloud is given in Fig. 28. Inside a cloud, the retarded cooling wave, which is 
associated with a high peak at the curve of cooling in the area of hydrogen 
helium recombination, is formed. An inverse temperature dependence of 
the cooling rate results in the development of a thermal instability and the 
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Figure 27: 27. Light curve of the optical afterglow in stellar magnitudes. 
The solid line gives the upper limit, and the dashed line the lower limit. The 
origin of coordinates corresponds to the moment of the GRB flare with an 
energy E and a total fluence near the Earth Fqrb — erg cm~^: curves 
la and lb correspond to the cases oi E — 10^^ erg, no = 10^ and E — 10^^ 
erg, no = 10^ respectively, [25]. 
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Figure 28: Temperature distribution in a cold homogeneous molecular cloud 

with the radius R = 1.5 pc, the hydrogen concentration uh = 10^ cm~^, 
the normal chemical composition after a homogeneous GRB flare of the total 
energy Q = 5 ■ 10^^ erg, with the flat spectrum for Emax = 2 MeV at the 
consecutive moments of time after the GRB: (1) 2.02, (2) 2.91, (3) 3.14, (4) 
3.70, and (5) 4.26 yr [10]. 

appearance of a thick cold layer between the hot center and thin external 
heated layer. It is interesting to note that the phase velocity of the external 
boundary of the cooling wave exceeds the velocity of light, which leads to 
the con- flnement of photons produced in the recombination inside the cold 
layer [25, 10]. 

The molecular cloud can be nonuniform, and the GRB flare can take 
place generally outside its center. It is conventionally assumed that the 
GRB radiation should be anisotropic being a narrow beam. The results of 
calculations of the distribution of density for such a general case are shown 
in Fig. 29. Here, the cooling proceeds faster in the higher-density region 
so that a thin hot layer is formed behind the GRB-propagation front, and a 
much thicker hot layer remains in the lower-density region behind the GRB. 
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Figure 29: Temperature distribution (dark filling, logT is indicated in the 
column on the right) and density levels (lines) after the GRB in the cold 
molecular cloud of a radius R— 1.5 pc with a hydrogen concentration uh — 

_ 2 / 2 

lO^e ^ /'"ocm~^, ro = 0, 75 pc, with the normal chemical composition after 
an anisotropic GRB with the angular distribution of radiated energy in the 

-(-)' 

form Q — 5- lO^^e V^o/ ^ _ q 4 j-g^^j^ g^^^j ^ f^^^ spectrum with £^^00; = 2 
MeV located at the distance xo — 1,0 pc from the molecular-cloud center. 
All distances are indicated in parsecs and the time in years [10]. 
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8.2 Hydrodynamic motions of gas caused by cosmo- 
logical GRB 

The maximal velocities attained by the matter of a homogeneous molecular 
cloud after the GRB flare inside it with a total energy release of 1.6 • 10^"^ or 
1.6-10^^ erg and the maximal photon energy E^ax — 2 MeV in a flat spectrum 
are equal to 2, 16 • 10^ and 5, 17 • 10^ km/s, respectively. The further increase 
in the GRB energy does not lead to increasing the velocity behind the shock- 
wave front because of the Compton-heating properties [10]. Fast cooling of 
matter hampers Shockwave enhancement even in the case of a nonuniform 
cloud with a moderate density gradient. 

It is possible to achieve a higher acceleration of matter after the GRB flare 
in some cases of a special distribution of matter around the GRB. When the 
GRB flares inside a low-density cone inside the molecular cloud, the density 
gradient is formed inside this cone and results in the jet of matter with a 
high velocity along the cone axis. The molecular-cloud matter nearby the 
GRB including that in the cone is heated to almost identical temperatures 
so that the arising pressure gradient follows the initial density gradient. The 
motion of matter to the cone axis under the action of this gradient results 
in the appearance of a cumulative jet along the cone axis, the matter in 
which is accelerated to the velocity of 1, 7 • 10^ km/s. The evolution of the 
temperature distribution and the induced velocity fleld are shown in Figs. 
30 and 31. The conic-type density distribution arises when the GRB flares 
between molecular clouds in the low-density region (in this case, the jet is 
ringshaped), or the similar cone can be formed by an anisotropic star wind 
of a GRB-progenitor star. 

The afterglows from an anisotropic GRB flare can have a strong angular 
dependence because of a strong absorption of light by dust, which is not 
evaporated after the flare along the directions perpendicular to the GRB 
axis. This fact can weaken the restrictions on a coUimation angle following 
from the lack of optical orphan bursts because the optical afterglow of similar 
type is coUimated almost the same as the GRB itself. The reradiation on the 
dust leads to the appearance of an infrared GRB afterglow, which should be 
more isotropic than the optical one. Therefore, the search for infrared orphan 
bursts accompanying an anisotropic GRB seems to be more promising. 

9 Afterglows in high-energy region 

Observations in the EGRET experiment onboard of the orbital Compton -ray 
observatory (CGRO) show that the GRBs radiate also in hard gamma regions 
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Figure 30: Distributions of velocity (arrows) and temperature (dark filling, 
logT arc indicated in the column on the right) in the initially cold molecular 
cloud of the radius i? = 1.5 pc at various moments of time after an isotropic 
GRB for the variant with an initial hydrogen- concentration distribution of 
the conic type uh = I0^e~^~^™***^^°^) cm~^. The calculated region is limited 
to the angles < 9 < n/lO, while vg = 0, on the boundary 6 = vr/lO; the 
distance from a GRB source is indicated in pc, and the timein yr [10]. 
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Figure 31: Zoomed central part of Fig. 30 [10]. 



up to an energy of 20 GeV [35]. Hard 7 rays were detected approximately for 
10 GRBs, photons with energies over 100 MeV being observed for 5 GRBs 
listed in Table 3 from [93]. As a rule, the hard 7 rays last longer than the 
basic soft GRB by up to 1.5 h in GRB 940217 (see Fig. 32). The comparison 
of the angular aperture in the experiments EGRET and BATSE and also 
the duration of their operation at the orbit results in the conclusion that 
the hard 7 rays should be observed during a significant part (approximately 
one third to half) of the entire GRB. The slope of the spectrum in the hard 
gamma region is in the range between 2 and 3.7 and quickly varies becoming 
softer with time (GRB 920622 from [94]). According to the observations in 
the hard gamma region of a radio pulsar in the Crab nebula [73] and PSR 
B1055-52 [102], the slope and the variation of their spectra have properties 
similar to those of a GRB. Taking into account the nonpulsing component of 
the Crab-nebula radiation, the slope of its gamma spectrum varies between 
1.78 and 2.75. If a GRB is associated with the supernova explosion and the 
neutron-star birth, the residual neutron-star vibrations after explosion can 
be responsible for a long afterglow in the hard gamma region [16, 103, 33]. 

All existing projects in the hard gamma region assume the investigation 
of GRBs as the primary problem. For this purpose, it is required to have the 
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Table 3: Observations of powerful GRBs in the experiment EGRET from 
[93] 



GRB name 


Maximal energy 
(GeV) 


Radiation 
duration 


Spectral 
function 


Delayed 
radiation 


GRB910503 


10 


84 s 


^-2.2 


X 


GRB910601 


0.314 


200 s 




X 


GRB930131 


1.2 


100 s 


^-2.0 


X 


GRB940217 


18 


1.5 h 


^-2.6 


X 


GRB940301 


0.16 


30 s 


^-2.5 






Figure 32: GRB flare of February 17, 1994, the gamma rays from which 
were observed in GeV-energy region in the experiment EGRET within 1.5 
h after the basic flare. The composite figure includes the EGRET data and 
the data of observations of the basic flare in the BATSE experiment and at 
the sateUite Ulysses [35] . 
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angular aperture as large as possible. It is expected that the nearest-in-time 
launch will be that of the small astronomical gamma sateUite of the Itahan 
Cosmic Agency (AGILE) of about 130 kg in weight, which should operate 
within 3 years in the energy range from 30 MeV to 30 GeV according to the 
project. In the AGILE experiment, it is proposed to use solid-state silicon 
detectors the fabrication technology of which is well developed. It is proposed 
to achieve the optimum angular resolution when constructing the images in 
the gamma region at a level of 5' - 20' for powerful sources and to perform the 
observation in an unprecedented large field of view of about 3 sr. This results 
in an effective sensitivity which is comparable with that of the EGRET for 
point sources observed along its axis and much higher for the sources away 
from the axis. The detection of GRB is expected in the range above 50 MeV 
at a level of 5 - 10 events per year [1]. 

The gamma large-area space telescope (GLAST), slated for launch later, 
is an international space project with the participation of the space agencies 
of various countries. The project is purposed for the space investigations in 
the energy range from 10 keV to 300 GeV. After the reahzation of several 
successful research space projects in the field of gamma-ray astronomy, the 
EGRET project was realized onboard of the CGRO launched in 1991. In 
the EGRET experiment, the first complete survey of the sky in the energy 
range from 30 MeV to 10 GeV was performed, and the majority of sources 
discovered in this survey remained unidentified. In the GLAST experiment, 
the investigations of the gamma region will be carried out in the range 20 
MeV - 300 GeV. The comparative characteristics of GLAST (EGRET) are as 
follows. The maximal effective region amounts to 8000 cm^ (1500 cm^), the 
field of view exceeds 2 sr (0.5 sr), the angular resolution is better than 3.5° at 
100 MeV and better than 0.15° at energies higher than 10 GeV (5.8° at 100 
MeV), the energy resolution is better than 10% (10%), and the localization 
accuracy (radius at the l-a level) for steady high-altitude sources with large 
\b\ and a flux of 10^'' cm^^s^^ in the energy range higher than 100 MeV 
amounts to less than 0.5' (15'). On board GLAST, the almost-all-sky monitor 
in the soft gamma region is specially mounted to obtain simultaneously the 
data in soft and hard gamma-photon energy regions for the independent 
detection of GRBs. As in the experiment AGILE, the detector for detecting 
hard gamma photons comprises 18 silicon strips. This technology has a long 
and successful history of application in experimental high-energy physics in 
accelerators. It has a high detection efficiency (>99%), an excellent positional 
resolution (<60 //m in this design), and a high signal-to- noise ratio (>20 : 
1). GLAST is expected to have high detection efficiency at energies above 
10 GeV and the possibility of GRB localization with accuracy reasonably 
high for the fast search of their radiation at all larger wavelengths. GLAST 
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is expected to detect about 200 GRBs per year, half of them located in the 
radius of no more than 10 arcmin, giving a good possibility to obtain the 
image of the region with GRBs using a large-aperture optical telescope [72] . 

The original design of the gamma telescope in the energy region higher 
than 0.5 - 1 GeV was considered in [24, 58, 31] and named CYGAMCylindri- 
cal Gamma Monitor. The principal schematic drawing of GYGAM is shown 
in Fig. 33. The hollow inside cylindrical shape allows the authors to reduce 
the weight of the device several times for the same sensitivity and the best 
angular resolution at an energy of about 1 GeV. In this scheme, the drift 
chambers [31] were stipulated for detectors; however, the above silicon strips 
can also be used instead of them. It is possible to expect that, in the energy 
range of 1 - 10 GeV, the GYGAM should have characteristics comparable 
to or even exceeding those of the GLAST, being almost three times lighter 
and, correspondingly, cheaper to manufacture. The large field of view of 
approximately 7 sr will be highly suitable for detecting GRBs. 

10 Lines in the hard x-ray region 

Lines in the hard x-ray region of GRB spectra were discovered in the experi- 
ment KONUS [64] . They were interpreted as cyclotron lines and observed in 
approximately 30% of GRBs. The lines and spectra were strongly variable: 
an observed absorption-line depth decreased in time (Fig. 34). The BATSE 
detector had the lower spectral resolution, and, for a very long time, no lines 
were found in the observed spectra. In later publications, at last, the pres- 
ence of spectral details in the hard-x-ray region of BATSE GRB spectra was 
reported [29]. In this paper, when investigating the spectra of 117 GRBs, 
statistically significant lines were found in 13 GRB spectra. A good example 
is GRB941017 for which the consistent lines are well seen in the spectra of 
two detectors. For some GRBs, one detector observed a line, while it proved 
to be statistically insignificant for another detector. Figure 35 shows the 
spectrum of GRB 930916 with the statistically significant line from [29]. In 
this paper, the authors affirmed that the doubts about the reality of the 
spectral details of GRBs in the hard x-ray region will remain until the causes 
for the differences between spectra of the same GRB observed by different 
detectors are understood. However, we note that the GRB 930916 spectra 
were obtained only within 20 s after the flare onset, and, according to [64], 
the strongest spectral details were observed at the GRB initial stage (see 
Fig. 34). The only study devoted to the interpretation of such lines in the 
cosmological model [40] is based on considering the spectrum strongly shifted 
in the blue region (the corresponding relativistic factor F = 25 — 100) for the 
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Figure 33: Schematic diagram of proposed CYlindrical GAmma Monitor 
(CYGAM). ACi are the coincidenceanticoincidence counters, ABi is the butt- 
end anticoincidence shield, W is the lumped converter, Ci are the coincidence 
counters; and DCi are the drift-chamber layers [24]. 
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Figure 34: Spectral evolution of GRB of November 1, 1979. (12) The spec- 
trum obtained during the first 8 s with the absorption line at an energy of 
~ 65 keV and a wide emission structure in the region of 350 - 650 keV, (4) 
the spectrum of the fourth 4-s interval, and (68) the total spectrum of the 
sixth, seventh, and eighth intervals [65]. 

gas cloud illuminated by the fireball 7 rays. A similar model was proposed 
in [21] for an explanation of the lines observed in the experiment KONUS 
within the framework of the GRB galactic model [23]. The model was based 
on the explosion near the neutron-star surface resulting in the jet of an ex- 
panding cloud with v/c = 0.1 — 0.3 in which a line absorption occurs (see 
Fig. 36). In this model, the observed attenuation in the absorption versus 
time (Fig. 34) was related to a decrease of the absorption thickness during 
the expansion of the cloud. 
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Figure 35: Spectrum of GRB 930916 obtained from the BATSE observations 
in the interval from 22.144 to 83.200 s after its detection. The bump at 30 
keV is associated with the K-boundary of iodine in the Nal detector. On the 
left: the best fit of data from the detector SD 2 with the continuum spectrum. 
On the right: with the addition of a narrow spectral line of radiation at 45 
keV, the value is improved by 23.1. Solid line corresponds to the total 
spectrum, and strokes separately give the continuum spectrum and the line 
[29]. 
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Figure 36: Time evolution of the absorption line formed due to the presence 
of barium ions with Xsa = 1/300 in the iron plasma of an expanding cloud 
located between a GRB source and an observer, ta < h < tc < td, is the 
GRB spectrum before the passage through a cloud [21]. 
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11 GRB and origin of ultrahigh-energy cos- 
mic rays (UHECRS) 

In [69], it was noted that the squares of errors of two events associated 
with cosmic rays (CRs) of the highest energies overlap with the squares of 
errors of strong GRBs. GRB 910503 and GRB 920617 were observed 5.5 
and 11 months before the corresponding CRs, respectively. In one of the 
cases, the most powerful CRs coincide with the most powerful GRB from 
the BATSE catalogue. Accidental coincidences of such events is very low; 
therefore, it was assumed in [69, 70] that the same phenomenon results in the 
occurrence of both GRBs and ultrahigh-energy CRs. A time delay, as well as 
a small divergence at the celestial sphere is expected because the magnetic 
field affects the trajectory of a CR charged particle and makes it curvilinear. 
The estimates show that only the Galaxy field taken alone can lead to an 
observed delay of CRs in comparison with GRBs. If we assume that the 
energy release in CRs and GRBs is approximately identical [69], several tens 
of percents of the strongest GRBs can explain the available experiments on 
detecting ultrahigh-energy CRs (UHECRs). 

The assumption that GRBs are the sources of UHECRs was almost si- 
multaneously made in [110, 111, 108]. The isotropy of distribution of both 
events in the sky and their cosmological origin point to the GRB UHECR 
relation. In addition to this, the observed energy-release rate (per unit vol- 
ume) in GRB 7 rays proves to be comparable with the high-energy-proton 
energy, which is necessary for explaining the observable UHECR fiux. The 
energy-release rate was estimated in [112] (see also [109]). It was assumed 
that the extragalactic protons in the energy range from 10^^ to 10^^ eV are 
produced in cosmologically distributed sources with the following density and 
spectrum: 

^ 0.6 X lO^^ergMpc-V^'- (H) 

dEp 

The energy-release rate in a GRB can be estimated assuming that a local {z = 
0) GRB frequency amounts to ~ 0.5Gpc~^yr~^ [92]. In this case, an average 
GRB energy with known redshifts in the MeV range of the 7-photon energy 
is equal to fa 0.5/ (Gpc^yr) erg. The GRB frequency with known redshifts is 
less than the total observed GRB frequency by a factor of ^ 0.7 x 0.5 = 0.35. 
The factor 0.7 is associated with the GRB detection threshold by the satellite 
Beppo-SAX, which is approximately two times higher than that found from 
the BATSE data used for estimating the total GRB frequency. The factor 
0.5 is associated with the fact that optical afterglows are detected only for 
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approximately half of GRBs. Thus, the local {z = 0) energy release rate in 
the MeV-photon region during the GRB flares is equal to 

S[MeV] ^ 0-35 X 3 X lO^^erg x 0.5/GpcVyr = 0.5 x lO^^erg/MpcVyr. (12) 

This energy-release rate is comparable with that for UHECRs from Eq. (11). 

The cncrgy-rclcasc rate can be estimated in another way using the lumi- 
nosity function obtained in [92] . On the basis of the large set of GRB events, 
it was found that the number of GRBs as a function of the peak luminosity 
Lc follows a power law in the range from ~ IQ^^-^ to ~ 10^^'^ crg/s and then 
abruptly falls. Using a value of ~ 10^^ erg/s for the peak luminosity and 10 
s as an average equivalent time of the total energy release, we estimate the 
total GRB energy-release rate to be approximately an order of magnitude 
less than that followed from Eq. (12). This value is more than ten times 
less than the UHECR energy-release rate from Eq. (11). With taking into 
account the large number of uncertainties, even the divergence in order of 
magnitude cannot be the basis for the final refusal to consider the hypothesis 
about the GRB - UHECR relation. 

12 Short GRBs and SGRs 

The cosmological GRB origin creates many problems when constructing a re- 
alistic physical model. The basic problem consists in the difficulty of obtain- 
ing an enormous energy release over such a short time in a small volume. For 
a coUimated GRB, the requirements on the energy release decrease; however, 
the coUimation angle cannot be too small. Observations of optical orphan 
bursts during the complete sky survey by optics can essentially improve the 
restrictions on the coUimation angle. The measurements of optical spectra 
of prompt afterglows, while optical luminosity is still reasonably high, and 
the investigation of the polarization in optical and x-ray regions are very 
important at this stage for clarifying the radiation mechanism. 

It is probable that GRBs do not represent a homogeneous set of events 
and consist in bursts of completely different origins. The statistical analy- 
sis reveals at least two separate sets consisting of long-duration (>2 s) and 
short-duration bursts. We note that optical afterglows with determining the 
redshifts were measured only for long bursts. Therefore, it is possible that 
the short bursts are of another origin (possibly, even galactic). 

It is interesting to compare the properties of short GRBs with giant flares 
of soft gamma repeaters (SGRs) located inside the Galaxy. If the SGRs were 
located at a greater distance on which a usual activity of SGRs would not be 
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seen, and only giant bursts would be detected, these last would be attributed, 
without doubt, to conventional short GRBs. The hght curves of two giant 
bursts are shown in Figs. 37 and 38 according to the observations by the 
KONUS-WIND device [66, 67]. The frequency and power of giant SGR bursts 
observed already in all four reliably identified SGRs in the Galaxy and in the 
Large Magellanic Cloud (LMC) are such that there is a real possibility of 
their observations from local-group adjacent galaxies as short GRBs. The 
estimates show that, by statistics, more than 10 short GRBs of similar type 
should be observed from the Andromeda Nebula (M 31) and other local-group 
galaxies [18]. No GRBs in nearby adjacent galaxies could point to the fact 
that the observed SGRs are closer and weaker objects than it is believed now 
on the basis of their possible genetic relation to supernova remnants (SRs), 
the distances to which can be estimated. It is possible that this relation does 
not exist at all [18]. 

13 Conclusions 

Despite more than 30 years of observations and numerous attempts at con- 
structing theoretical models, until now, there is no authentic model either 
for the formation of observable radiation or for the GRB energy source. Al- 
though long GRBs are most likely of the cosmological origin, it is possible 
that short GRBs have another origin (probably, even galactic) and are asso- 
ciated with giant SGR flares [65] . 

Critical experiments, which could finally clarify the GRB nature, are as- 
sociated with observations in all regions of the electromagnetic spectrum. In 
optics (from infrared to ultraviolet), it is necessary to obtain GRB spectra 
at the earliest stages as close as possible to the time of the -ray burst it- 
self. This problem is formulated for wide-angle automatic telescopes, which 
already operate or are being designed in various countries. They should be 
quickly directed to a sky region using the information obtained from a gamma 
satellite right after the detection of a GRB. Similar problems is solved by the 
specialized satellite SWIFT launched on November 20, 2004 and purposed 
for observations of GRBs in a wide range from hard x rays (150 keV) to the 
optical band. SWIFT registers about 100 GRBs with a localization accu- 
racy of about 1' within one year of operation. An important observational 
problem in optics is the search for the orphan bursts not accompanied by a 
7-ray pulse, about the degree of collimation of the gamma pulse. Such obser- 
vations require a long-term all-sky scanning using optics; they arc possible 
with either a system of continuously operating very wide-angle telescopes 
covering the entire sky or a specialized satellite with an optical telescope 
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Figure 37: Light curve of the giant burst SGR 1900 + 14 of August 27, 1998. 
The presented data on the radiation in the photon-energy region £^ > 15 keV 
are from [66]. 
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Figure 38: Light curve of the giant burst SGR 1627-41 of June 18, 1998. The 
presented data correspond to the radiation in photon energy region E > 15 
keV. Time of rise of brightness amounts approximately to 100 ms [67]. 
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(or system of telescopes)the all-sky monitor. In addition to the gamma ob- 
servations for detecting GRBs in the energy region of about one MeV, the 
essential information can be obtained from observations in the hard gamma 
region with a photon energy ranging from a hundred MeV to tens of GeV. 
Such observations in the same energy range are planned to be carried out by 
the gamma large-angle space telescope (GLAST) with the expected launch 
in 2007 and by the gamma small space telescope AGILE, which should be 
launched earlier. The collapse of a high-mass star nucleus with the formation 
of a black hole surrounded by a massive rotating disk represents the probable 
cosmological GRB model. The fast fall of the disk onto the black hole as a re- 
sult of the magnetorotational processes accompanied by an enormous energy 
release in the form of neutrino, radiation, and electronpositron pairs can be 
the origin of GRBs. The possibility of originating GRBs in exotic processes 
associated with cosmic strings, antimatter, etc., is also not excluded. 
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